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Unipolar Brush Cells (UBCs) have been suggested to play a critical role in cerebellar 
functioning, yet the corresponding cellular mechanisms remain poorly understood. 
UBCs have recently been reported to generate, in addition to early-onset glutamate 
receptor-dependent synaptic responses, a late-onset response (LOR) composed of a slow 
depolarizing ramp followed by a spike burst (Locatelli et al., 2013). The LOR activates 
as a consequence of synaptic activity and involves an intracellular cascade modulating 
H- and TRP-current gating. In order to assess the LOR mechanisms, we have developed 
a UBC multi-compartmental model (including soma, dendrite, initial segment, and axon) 
incorporating biologically realistic representations of ionic currents and a cytoplasmic 
coupling mechanism regulating TRP and H channel gating. The model finely reproduced 
UBC responses to current injection, including a burst triggered by a low-threshold spike 
(LTS) sustained by CaLVA currents, a persistent discharge sustained by CaHVA currents, 
and a rebound burst following hyperpolarization sustained by H- and CaLVA-currents. 
Moreover, the model predicted that H- and TRP-current regulation was necessary and 
sufficient to generate the LOR and its dependence on the intensity and duration of mossy 
fiber activity. Therefore, the model showed that, using a basic set of ionic channels, UBCs 
generate a rich repertoire of bursts, which could effectively implement tunable delay-lines 
in the local microcircuit. 



Keywords: biorealistic modeling, unipolar brush cells, vestibular cerebellum, slow synaptic responses, ionic 
channel regulation 



INTRODUCTION 

The classical neurotransmission process is based on release of a 
chemical transmitter followed by activation of ionotropic recep- 
tors and the consequent generation of membrane potential 
changes in the postsynaptic neuron (Katz and Miledi, 1964). In 
addition to this, slow neurotransmission mechanisms have been 
reported, which involve the activation of intracellular second- 
messenger cascades and slow regulation of membrane ionic chan- 
nels. These can either be voltage-dependent channels taking part 
to regulation of neuronal electroresponsiveness, or be ionic cyclic 
nucleotide-gated channels (CNGCs) specifically involved in neu- 
romodulatory responses. While these latter are well known in 
heart cells and photoreceptors, they have only recently been rec- 
ognized in neurons (Brown et al, 1981; Wainger et al., 2001; 
Clapham, 2003). Interestingly, the unipolar brush cells (UBCs) of 
the cerebellum expresses both responses generated by ionotropic 
receptors (Rossi et al., 1995; van Dorp and De Zeeuw, 2014) 
and a late-onset response (LOR) generated by receptor coupling 
through cytoplasmic messengers to membrane ionic channels 
(Locatelh et aL, 2013). 

UBCs are interneurons of the cerebellar granular layer form- 
ing a complex network with mossy fibers, granule cells, Golgi 



cells, and other UBCs (Mugnaini et al, 1997). UBCs have been 
proposed to play an essential role for shifting and converting 
the phase of mossy fiber activity required for controlling and 
consolidating motor learning (Gao et al., 2012). When injected 
with depolarizing currents from a negative membrane potential, 
UBCs respond with a calcium spike priming a sodium spike burst 
(Russo et al, 2007). Then, as depolarizing currents are increased, 
the initial burst is followed by a tonic discharge (Diana et al., 
2007). The ionic currents responsible for UBC electrorespon- 
siveness have been partly clarified, including CaLVA and CaHVA 
currents, the H-current and different types of K currents (Diana 
et al., 2007). When activated synaptically, the LOR was shown to 
depend on the intracellular cAMP cascade and on subsequent reg- 
ulation of H- and TRP-currents (Locatelli et al., 2013). Despite 
these observations, it remained unclear whether the known set 
of UBC ionic currents was sufficient to explain membrane elec- 
trogenesis, including responses to current injection and the LOR, 
and whether an intracellular messenger could effectively couple 
the neurotransmission process to LOR generation. 

In this paper, intrinsic UBC electroresponsiveness was repro- 
duced in a biologically realistic UBC multi-compartmental 
model. The model could accurately predict all known aspects of 
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UBC electroresponiveness demonstrating the central role of H- 
and Ca^+ currents. The model then effectively predicted the LOR, 
when a generic intracellular mechanism was used to couple mossy 
fiber stimulation to TRP and H channel gating. LOR delay, dura- 
tion and spike frequency turned out to be efficiently regulated 
by the intensity and duration of synaptic stimulation. Therefore, 
the model confirmed that a fundamental set of ionic channels 
was sufficient to explain both UBC electroresponsiveness and 
the LOR. This mechanism could contribute to generate tunable 
delay-lines in the cerebellar network (Kennedy et al., 2014). 

METHODS 

This work presents a combined modeling (I) and experimental 
analysis (II) of UBC electroresponsiveness, in which the models 
was accurately matched to biological responses (III). The model 
was written in NEURON (Mines et al, 2007; Davison et al, 2009). 

THE UBC MODEL 

A multi-compartmental model representative of the main 
morphological components of a UBC was constructed using 
the NEURON simulator (NEURON version 7.3; Hines and 



Carnevale, 2001). The model consisted of soma, dendrite, initial 
segment and axon compartments generating a morpho-electrical 
equivalent of the UBC (Figure lA). The voltage- and Ca^+- 
dependent mechanisms (see below) were distributed among the 
compartments. With this approach, the model reproduced sat- 
isfactorily basic aspects of UBC electroresponsiveness elicited by 
somatic current injection. It should also be noted that, although 
the existence of UBC subtypes has been suggested based on his- 
tochemical analysis, the basic electroresponsive properties were 
homogeneous in a large majority of UBCs (Locatelli et al., 2013). 
Thus, we have reconstructed a canonical UBC model, which 
simulates the typical electrophysiological behavior of UBCs and 
lies within the scatter of physiological parameters values. The 
model was coupled to a recording electrode to reproduce real- 
istic current-clamp conditions and was calibrated on the cellular 
response corrected for —10 mV liquid-junction potential. 

The model ionic conductances were non-homogenously 
distributed over different compartments (Figure IB; Table 1). 
Gating kinetics were corrected using a Qio = 5 for IcaLVA 
activation (Destexhe et al., 1994) and a Qio = 3 for all 
the other currents according to the relation Q^g^™ '''^P'''^'' 




FIGURE 1 I UBC morphology and passive properties. (A) A UBC recorded 
in patch-clamp whole-cell configuration and observed in transmitted light 
microscopy. The cell contour has been identified and color-filled. The different 
parts of the cell have been identified and labeled. (B) Morphology and 
distribution of ionic channels in the UBC model. Same color code as in 
Figures 3-6, 8. (C) The trace shows a voltage-clamp recording of a current 
transient elicited in a UBC by a —10 mV voltage step (blue trace) and the 
corresponding model reconstruction using the same protocol (red trace). The 
transient is enlarged in the inset to show the close adherence of model to 



experimental responses. The plot shows the relationship between input 
resistance (Rjn, GQ) and input capacitance (Qn, pF) calculated from the fitting 
of current transients, both in experiments and in the model. Experimental 
points were taken from the UBC population recorded for this paper (n = 30) 
and from Locatelli et al. (2013); unpublished, n = 51 . The model prediction 
falls in the middle of the experimental data distribution. The values range 
reported by other works are shown as gray shadowed areas: a (Diana et al., 
2007), b (Russo et al., 2007, 2008), c (Locatelli et al., 2013), d (Birnstiel et al., 
2009). 
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(Gutfreund et al., 1995) to account for temperature differences 
between experimental recordings and the model (30°C). Nernst 
equilibrium potentials were pre-calculated from ionic concen- 
trations used in current-clamp recordings and maintained fixed, 
except for the Ca^+ equilibrium potential, which was updated 
during simulation according to the Goldman-Hodgkin-Katz 
(GHK) equation. The maximum ionic conductances were regu- 
lated to fit the UBC responses to various stimulations (Traub and 
Llinas, 1979; Traub et al, 1991; Vanier and Bower, 1999; Achard 
and De Schutter, 2006; Druckmann et al, 2007, 2008, 2011, 2013; 
Solinas et al., 2007a,b). Ionic channel gating was modeled follow- 
ing the Hodgkin and Huxley formulation or Markov- chains for 
multi-state transitions using mathematical methods reported pre- 
viously (D'Angelo et al, 2001; Nieus et al, 2006; Solinas et al, 
2007a,b) and are reported in Table 1. Membrane voltage was 
obtained as the time integral of the equation (Yamada et al., 1989). 



dV 
7t 



(1) 



Where V is membrane potential, C,„ membrane capacitance, gi 
are ionic conductances and V, reversal potentials (the subscript i 
indicates different channels), and i,„j is the injected current. The 
intracellular Ca^+ concentration, [Ca^+], was calculated through 
the equation. 



d[Ca^+] 



-ICa 



dt 



(IF * Ad) 



Pca ([Cfl2+] - [Cfl2+]o) (2) 



Where d is the depth of a shell adjacent to the cell surface of area 
tt) Pca determines the loss of Ca^+ ions from the shell approx- 
imating the effect of fluxes, ionic pumps, diffusion, and buffers 
(Traub and Llinas, 1979; McCormick and Huguenard, 1992; De 
Schutter and Smolen, 1998), and [Cfl^^"]o is resting [Ca^+]. 

Morphology 

Model morphology was reconstructed by keeping the dimensions 
of every compartment within the ranges reported in the literature. 
The model contains a spherical cell body (8 [xm diameter; Rossi 
et al, 1995; Morin et al, 2001; Kalinichenko and Okhotin, 2005; 
Russo et al., 2007), a short dendritic shaft (50 |im length and 
2 |xm diameter; Difio and Mugnaini, 2000) that originates from 
soma and connects to a comparatively large compartment named 
brush (25.25 (xm length and 10 [iva diameter; (Mugnaini et al., 
1997, 201 1) on its other end. The surface of the brush is extended 
by unfolding the three-dimensional structure of dendrites of real 
UBCs, which are enriched by many filopodia/dendrioles that 
form a structure reminiscent of a brush. The UBC axon was 
reproduced with five identical compartments (80 (im length and 
0.5 [LTCL diameter; Kalinichenko and Okhotin, 2005) that are seri- 
ally connected to the axon initial segment (AIS). The AIS was 
made of 5 compartments of equal length (0.5 |xm) and decreas- 
ing diameter moving toward the axon (from 3.2 to 0.8 |xm; 
Figure IB). The micropipette that was attached to the soma was 
simulated as a passive cable with sealed end, null capacitance, infi- 
nite wall resistance and an axial resistance (20 MQ,) equivalent to 
the sum of series and access resistance measured experimentally. 



Passive properties 

The axial specific resistance of axon and dendrites was set to 100 
Qcm, the specific membrane resistance was set to 47.6Kf2/cm^ 
and specific membrane capacitance was set at 1 |xF/cm^ (cf. 
Roth and Hausser, 2001; Diwakar et al, 2009, 2011). The input 
resistance (R;n, Gf2) and input capacitance (Cin, pF) of the 
model were calculated from the current transient elicited by 
250 ms voltage step of AV = — 10 mV from the holding potential 
of — 70mV (Figure IC). The charge displacement during volt- 
age steps followed a tri-phasic trajectory reflecting electro tonic 
compartmentalization, with the faster transient being related to 
the soma, the intermediate transient being mostly contributed 
by the brush and dendritic shaft and the slowest transient being 
mostly related to the axon (Locatelli et al., 2013); see Section II 
in Methods). Current transient analysis yielded Rin = 1.04 Gf2 
and Cin = 16.5 pF. The Cin value corresponds closely to the total 
cell capacitance of 16.7 pF calculated by summing the capaci- 
tance of the soma (2 pF), brush and dendritic shaft (12.7 pF), and 
axon (2 pF). 

Mechanisms of intrinsic electroresponsiveness 

The mechanisms of intrinsic electroresponsiveness were based 
on experimental observations including current-clamp (Locatelli 
et al., 2013) and voltage-clamp recordings (Diana et al., 2007; 
Birnstiel et al, 2009; Locatelli et al, 2013). Most currents intro- 
duced into the model were based on experimental observations 
like IcaLVA, IcaHVA, IkCu, Ih, Itrp, Inu while others includ- 
ing /fenj-, 7jfv, Ika, IksIow were introduced because needed to 
control the passive properties as well as the timing and fre- 
quency of spike discharge. In summary, the model included 10 
defined ionic mechanisms (see Table 1 for their mathematical 
representation): 

Ileah background leakage current. 
Itrp> cyclic nucleotide gated leakage current. 
Ih, hyperpolarization activated cationic current (H-type). 
IcaLVA, low-voltage activated calcium current (T-type). 
IcaHVA, high-voltage activated calcium current (L-type). 
lKCa> calcium-dependent K current current (BK-type). 
iMay Na current with transient, resurgent and persistent com- 
ponents. 

Ikv> delayed-rectifier K current. 

Ika, fast-activating K current (A-type). 

IksIow, slow- activating K current (M-type). 

Leakage currents (Iieak> Itrp)- The model included a leakage cur- 
rent (Iieak) adapted from a granule cell model (D'Angelo et al., 
2001) and localized in all 13 compartments of the UBC model 
with a reversal potential at — 70mV. The corresponding conduc- 
tance was tuned to obtain appropriate passive current transients. 
The presence of Ijnp in UBCs was reported by Russo et al. (2007) 
and confirmed by Locatelli et al. (2013). TRP channels in UBCs 
are ruthenium red sensitive, voltage-independent cationic cur- 
rents. TRP channels were modeled as a simple leakage channels 
with a reversal potential at OmV (Petersson et al., 2011). The 
model was adapted to change its conductance in response to 
second messenger concentration (see below). 
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H current (Ih)- The presence of Ih in UBC was reported by 
Russo et al. (2007) and confirmed by Locatelli et al. (2013). Ih 
in UBCs is manifest during the application of hyperpolarizing 
current steps causing sagging inward rectification (Santoro et al., 
2000). 7h proved able to slowly depolarize the UBC over the whole 
subthreshold range and causing a steady depolarization (~4 mV) 
between —50 and — 70mV (Locatelli et al, 2013). After regulat- 
ing current densities to match the sag and the steady-state level of 
hyperpolarizing responses, the model generated a robust rebound 
depolarization, as observed experimentally. The Ih model was 
adapted from the HCN current in thalamo-cortical neurons 
(McCormick and Pape, 1990) and was tuned using experimen- 
tal data on maximum conductance and channel kinetics in UBCs 
(Locatelli et al, 2013). The reconstructed Ih showed a 10% acti- 
vation at — 108 mV and half-activation at —91.5 mV with a slope 
factor of 7.7 mV. The Ih activation time constant was mono- 
exponential and decreased toward negative membrane potentials 
with mean values ranging from ~ 180 ms at — 130 mV to ~630 ms 
at — 90mV, in keeping with experimental results (Locatelli et al., 
2013). The model was adapted to shift its activation curve in 
reaction to second messenger concentration (see below). 

Na currents (Inu)- UBCs express TTX sensitive voltage-gated 
channels with transient, persistent and resurgent current com- 
ponents (Afshari et al., 2004). The Na channels of these neurons 
tend to recover from inactivation through open states, giving rise 
to a "resurgent" current that flows upon repolarization from pos- 
itive potentials (Raman and Bean, 1997; Afshari et al., 2004). The 
recovery from inactivation is associated with the resurgent cur- 
rent, which shortens the refractory period between action poten- 
tials favoring the generation of high-frequency bursts (Raman 
and Bean, 2001; Khaliq et al., 2003). As the resurgent properties 
are typical of Navl.6 (Raman and Bean, 1997; Afshari et al, 2004), 
we adopted the corresponding granule cell model (Magistretti 
et al., 2006). These channels have been inserted in soma, hillock 
and axonal compartments according to the most common local- 
ization found in other neurons, e.g., granule cells (Diwakar et al., 
2009) and Purkinje cells (Palmer et al, 2010). 

K currents (Ikv> Ika> IksIow Ikco)- Multiple classes of potassium 
channels (K+) are inserted in the model. delayed rectifier 
(KV), the fast inactivating A-type (KA), the persistent M-type 
(Kslow), along with the Ca^+ -dependent K channels (BK-type) 
were adapted without any modifications from the granule cell 
model (D'Angelo et al., 2001; Nieus et al, 2006). A-type and 
BK-type calcium dependent potassium channels were localized 
in brush, dendrites (McKay et al., 2005) and soma compart- 
ments with different maximum conductances to regulate first 
spike delay, spike after-hyperpolarization (AHP) and interspike 
interval. KV was localized in soma, hillock and axonal segments 
to regulate spike AHP. A slow repolarizing M-type (IksIow) and 
a fast-activating A-type {Ika) K currents in UBCs cells have not 
been reported yet, but were needed to regulate spike delay and 
tonic discharge frequency. 

Cfl currents (IcaLVA> IcaHVA)- In UBC two types of calcium chan- 
nels are expressed and play an important role in controlling 



electroresponsiveness. UBCs are bimodal, as they can either 
fire high-frequency spike burst when stimulated from hyperpo- 
larized potentials or generate tonic discharge during sustained 
depolarization. Fast inactivating Ca^+ channels (CaLVA or T- 
type) generate low-threshold spikes (LTS) triggering the high- 
frequency spike bursts, while poorly-inactivating Ca^+ channels 
(CaHVA or L-type) sustain the tonic firing. T-type calcium 
channels in UBCs are of the Ca3.1 subtype (Llildebrand et al., 
2009). Channel kinetics were adapted from a previous model 
(Anwar et al., 2012) in order to reconstruct UBC specific cur- 
rent kinetics, voltage-dependence and time-dependence (Diana 
et al, 2007; Birnstiel et al, 2009). L-type calcium currents were 
also reconstructed adapting a previous model (D'Angelo et al., 
2001) accounting for UBC specific properties (Diana et al, 
2007; Birnstiel et al., 2009). In the case of T-type channels, the 
ionic conductance was computed using GHK equation and the 
Ca^+ permeability (Hirschberg et al, 1998; Anwar et al, 2012). 
Both types of calcium channels were localized in brush, den- 
dritic shaft and soma compartments. Maximum conductance of 
T-type channels was set at a higher level in brush compared 
to soma, whereas that of L-type channels was maintained the 
same in all compartments (Diana et al., 2007; Birnstiel et al., 
2009). 

Ca^'^ dynamics. Ca^+ dynamics were adapted to yield Ca^+ 
transients of about 1 (jlM following spike-dependent opening of 
CaHVA channels (Gabbiani et al, 1994). Parameters used in 
Equation (6) (d = 200 nm, Pca = 1-5, and [Ca^+Jo = 100 nM) 
were similar to those used in other neurons (e.g., Traub and 
Llinas, 1979; Traub, 1982; De Schutter and Smolen, 1998; 
D'Angelo et al, 2001). The computed intracellular Ca^+ concen- 
tration [Cfl^+]i activated I^ca just after the spike upstroke and 
generated the fast AHP. It should be noted that, in the model, 
IcaLVA was not allowed to contribute to the Ca^+ pool controlling 
Ikco- 

Special issues in modeling UBC electroresponsiveness. In UBCs, 
Ih, IcaLVA, and IcaHVA provided strong constraints for model 
tuning, since their conductances were known from experimental 
measurements, (i) After setting morphology and passive proper- 
ties, Ih was added with known conductance and its effect on the 
hyperpolarizing responses was evaluated, (ii) Then, IcaLVA and 
IcaHVA were added respecting the conductance, localization and 
proportions determined experimentally, (iii) The complement of 
spike-related currents was reconstructed with mechanisms taken 
from a previous granule cell model (D'Angelo et al., 2001; Nieus 
et al., 2006). The density of the Na+ current was adjusted to 
reproduce the action potential threshold and overshoot. The 
delayed rectifier K+ current, Ikv, was needed to generate spike 
repolarization and balanced with I^ca to match the spike shape. 
The requirement of a Ca-^+ -dependent voltage-dependent K+ 
current, Ikcu, emerged from the inability of Ikv to fully account 
for the fast phase of spike AHP and to guarantee Na+ channel re- 
priming at high discharge frequency. Ika was needed to regulate 
spike delay in response to depolarizing current injection and dur- 
ing rebound excitation and IksIow was required to regulate tonic 
discharge frequency. 
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Receptor coupling to ionic channels: modeling the late-onset 
response 

The LOR was reproduced in the model by modulating the In and 
IjRP ionic currents with respect to the amount of intracellular 
factors related to the cAMP cascade (Locatelli et al, 2013). As a 
matter of fact, the exact intracellular-biochemical cascade driv- 
ing the LOR is unknown. Therefore, we implemented a generic 
second messenger pathway reflecting general principles of bio- 
chemical organization (Wainger et al., 2001; Zhou et al, 2009; 
Shen et al., 2011). This includes a receptor X placed in prox- 
imity of the UBC synapse, which activates into X* over a slow 
time course (>200 ms; consistent with Locatelli et al., 2013) upon 
mossy fiber bundle stimulation (Equation 3). Then X* activates a 
second-messenger cascade transforming Y into Y* (Equation 4). 
Finally, Y* modulates the shift of Ih activation curve and the TRP 
channel gating depolarizing the model and triggering the LOR. 
The cascade reaction is the following (see Figure 4): 

a 

X±^X* (3) 
S(X*) 

Y ^ Y* (4) 
Y 

Where, 

nspk 

a = ^ e (f - 0 © (ts + r-t)co (5) 

s= 1 

S (X*) = l/(exp ((X* - X*_half)/K)) (6) 

Where, a and f} are the rates of forward and backward reactions, 
a is a step function that, at the time of occurrence (tj) of each 
input pulse, increases by the rate of synaptic weight (w) and then 
decays back after a time interval r (Equation 5). yS is a backward 
rate constant that determines a fixed time to degrade X* into X. 
The initial value of X was calibrated in order to reproduce the 
LOR generated by a 10-pulse burst at 100 Hz. Formation of Y* 
from y is a reversible reaction (Equation 6). The forward reac- 
tion increases Y* depending on the amount X*. The initial Y 
value used in the model (5 mM) is in the range of the expected 
ATP level in neurons (Fitz, 2007). S and / determine the forward 
(Y Y*) and backward {Y* Y) reaction rates, respectively. 5 
is a sigmoidal function, with half-activation X*_half = 2.7 mM 
and slope K = 0.75 mM"^ 

The particle Y* was used to modulate the H and TRP channels. 
In the case of TRP channels, the Gmax-TRP value was increased 
proportionally to Y* (see Figure 4) (Equation 7). This mecha- 
nism is supported by observations of TRP channel modulation by 
intracellular cAMP (Petersson et al., 2011) and is consistent with 
the effects observed in UBCs (Locatelli et al, 2013). In the case 
of H channels, Y* modified In gating by proportionately shifting 
the steady-state activation curves to more positive potentials. This 
effect, observed in other cell types (DiFrancesco and Mangoni, 



1994), is consistent with that observed in UBCs (see Figure 4; 
Locatelli et al., 2013). Thus, the H channel half-activation value 
was shifted proportionately to Y* between —91.5 and —66.5 mV 
(Equation 8). 

new Gmax-TRP = Gmax-TRP * (0.16 -|- Y*) (7) 

newVi„act-H = Vinact-H + w( — jz-^ ^TSTTTTTt) 

Vl -I- exp(0.6 - Y*/0.i)/ 

WHOLE-CELL RECORDINGS FROM UBCs IN ACUTE CEREBELLAR SLICES 

UBC recordings were carried out as reported previously (Locatelli 
et al., 2013). Whole-cell patch-clamp recordings were performed 
from 30 UBCs in the internal granular layer of lobule X in acute 
rat cerebellar slices of P18-P25 wistar rats. Out of these, 12 UBCs 
were used to study intrinsic electro responsiveness and 18 to study 
the LOR. In addition, to estimate the passive properties and the 
resting membrane potential, data were also re-analyzed from a 
population of 51 UBCs reported by Locatelli et al. (2013) for a 
total of 81 measurements. 

In brief the patch pipettes were pulled with a horizontal 
puller (Sutter Instruments, Novate, CA, USA) from thick-walled 
borosilicate glass capillaries (Hingelberg, Germany) and had a 
resistance of 7-10 MQ, when filled with the intracellular solution 
(in mM): 126 potassium gluconate, 4 NaCl, 15 glucose, 5 Hepes, 
1 MgS04.7H20, 0.1 BAPTA, 3 ATP, 100 [xm GTP; pH adjusted to 
7.2 with KOH (in solution at 0.1 mM). The BAPTA-Ca2+ buffer 
was prepared as explained previously (D'Angelo et al, 1995; Gall 
et al, 2003). To minimize pipette tip capacitance, pipettes were 
coated with Sylgard (Dow Chemical, Midland, MI, USA) and 
the bath level was kept as low as possible. After careful pipette 
capacitance cancelation (D'Angelo et al., 1995), — lOmV 250 ms 
voltage step were applied to the cell from the holding poten- 
tial of — 70mV and the corresponding current transients were 
recorded (sampling frequency at 20-40 kHz, low-pass filtering 
5 kHz). The current transient decayed multi-exponentially, prob- 
ably reflecting charging of compartments corresponding to the 
UBC soma, dendritic shaft, brush, and axon. Tri-exponential fit- 
tings allowed us to estimate the electrode series resistance (Rs) the 
cell input capacitance (C,„) and the input resistance (Ri„) using 
the equation: 

7(f) =AiXexp(fc-f/Tl)-|-A2Xexp(fc-f/T2) (9) 
-FA3Xexp(/:-f/r3)-|-C 

Series resistance (Rs) was calculated as Rg = tvc/Qn, where 
Tvc = is the decay time constant of the current tran- 
sient related to somatic charging. In our extended UBC sam- 
ple, we measured Q„ = 16.9 ± 0.67 pF, _R,„ = 0.87 ± 0.07 Gf2, 
Rs = 22 ± 0.96 MQ (mean ± MSE, n = 81 for all parameters). 

All current-clamp recordings were performed using the fast 
current-clamp mode of the amplifier to accelerate membrane 
charging (D'Angelo et al, 1995, 1998; Prestori et al, 2008). The 
resting membrane potential was measured over 50 ms in 12 con- 
secutive traces. Intrinsic excitability was investigated by setting 
resting membrane potential at —80 mV and injecting 800 ms cur- 
rent steps (from —16 to 48 pA in 4-8 pA increments). Current 
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Steps were applied from the holding potential of — 80mV and 
the corresponding voltage deflections were recorded (sampling 
frequency at 20-40 kHz, low-pass filtering 5 kHz). 

Electrical stimulation was performed by placing a bipolar 
tungsten electrode over the mossy fiber bundle in the granular 
layer (D'Angelo et al., 1993, 1995). The stimuli consisted of volt- 
age pulses (0.25 ms, 5-15 V) organized in 100 Hz trains of varying 
intensity and length. 

DATA ANALYSIS AND MODEL MATCHING TO EXPERIMENTAL DATA 

Spike analysis 

The analysis of experimental and simulated traces was per- 
formed using identical measurement procedures and the pro- 
cess was automated using dedicated scripts written in MATLAB 
(MathWorks, Natick, MA, USA). In response to depolarizing 
currents, spike threshold was detected at the point where the 
depolarization rate reached 5 mV/ms (Forti et al., 2006). For each 
spike, the overshoot and its position were measured. The fre- 
quency ratio {instantaneous frequency I steady -state frequency) in 
response to depolarizing currents was calculated as the ratio 
between the first and fourth ISI, which are representative of the 
burst and tonic discharge regions, respectively. The sag ampli- 
tude was calculated as the difference between the minimum and 
steady-state value of membrane potential during responses to 
hyperpolarizing currents. When the LOR was elicited, its delay, 
duration and all the firing parameters were measured. Data are 
reported as mean ± s.e.m. 

Model construction, tuning, validation 

The modeling procedure could be divided into three phases: con- 
struction, tuning, and validation. In the UBC model, after defining 
the morphological and passive properties, setting the ionic chan- 
nel complement and the ionic channel gating properties (con- 
struction), the only free parameters remained the maximum ionic 
conductances. These were pre-adapted based on experimental 
estimates (see Figure 3) and then fine-tuned against experimen- 
tal data using voltage traces elicited in response to step-current 
injection (pulses from different holding potentials and responses 
to hyperpolarization) [tuning). Finally, the model was used to 
predict response properties, which were not included in the 
data-set used for construction (validation). Model validation was 
performed at several levels by evaluating the matching with: 

(i) Passive properties of a UBC population more extended than 
the one used to construct the model, 

(ii) Active electroresponsive properties comprising resting 
membrane potential, transition from burst to tonic firing, 
sagging responses to hyperpolarization, rebound excitation 
and burst inactivation from depolarized membrane poten- 
tials, 

(iii) LOR properties including dependence on the input pattern, 
intensity of the LOR generating currents, LOR pharmaco- 
logical sensitivity. 

Model robustness 

As quantitative validation, we assessed the ability of the model 
in maintaining typical UBC properties once one of the max- 
imum ionic conductances was varied in turn (Solinas et al.. 



2007a,b). Model robustness was evaluated both for eletrorespon- 
siveness (see Figure 6A), for the LOR (see Figure 6B), and for 
their combination (Figure 6C). 

RESULTS 

UBC MULTI-COMPARTMENTAL MODEL: CHANNEL DISTRIBUTION AND 
PASSIVE PROPERTIES 

UBCs (Figure lA) were recorded from the granular layer of lam- 
ina X of the rat vestibular cerebellum (Mugnaini and Floris, 
1994) using patch-clamp techniques in acute slices and were 
modeled using the NEURON simulator. In order to investigate 
whether the proposed LOR mechanisms were consistent with 
current knowledge on UBC intrinsic electroresponsiveness and 
whether the prediction of LOR generating mechanisms was cor- 
rect (Locatelli et al., 2013), we have developed a detailed UBC 
computational model, in which synaptic transmission was cou- 
pled to H- and TRP-channels through a cytoplasmic modulatory 
mechanism. 

The model was constructed respecting the morphological con- 
straints and passive properties (Figure 1) and using the available 
electrophysiological data on ionic channels of UBCs (Table 1). 
Then the model was calibrated through a multi-parametric 
comparison with current-clamp response to current injection 
(Figures 2, 3). Finally, the model was extended through a mech- 
anism coupling synaptic activity to ionic channel gating and 
became able to generate the LOR (Figure 4). The robustness of 
this procedure was evaluated by performing various tests allow- 
ing model predictions to be compared to native UBC responses 
(Figures 5-8). 

The UBC model (Figure IB) was composed of 4 main com- 
partments (brush, dendritic shaft, soma, axonal initial segment, 
axon). The model was endowed with 8 voltage-gated or Ca^^- 
gated ionic currents (Inu, hv, Ika, IksIow, ^KCa, IcaLVA, IcaHVA, W 
and with a voltage-independent leakage current (Itrp)- The kinet- 
ics of individual ionic currents were fitted to UBC data (when 
directly available) or adapted from previous cerebellar neuron 
models (D'Angelo et al, 2001; Anwar et al, 2012). The density 
of ionic channels in different model compartments was pre-set 
according to the known channel localization (e.g., Na channels 
were more expressed in the AIS, CaLVA channels were more abun- 
dant in the brush than in the soma while CaHVA channels were 
equally distributed in soma and brush; (Khaliq and Raman, 2006; 
Diana et al., 2007) and subsequent adjustments to the current 
density were made to properly fit the UBC response to current 
injection (cf. Solinas et al., 2007a,b; see below). 

The calibration of UBC passive properties was done with 
respect to electrophysiological measurements (Figure IC). After 
having reconstructed the UBC compartmental morphology, the 
injection of voltage step into the model yielded input capacitance 
and resistance values falling in the middle of the experimental 
data distribution (the average experimental values were 16.3 pF 
and 0.97 GQ, the average model values were 16.7 pF and 1 Gf2). 

UBC ELECTRORESPONSIVENESS AND THE LATE-ONSET RESPONSE 
(LOR) 

In current-clamp whole-cell recordings, UBCs showed a typi- 
cal electroresponsive pattern comprising the following elements 
(Figure lA): 
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FIGURE 2 I Recordings and simulations of UBC responses. The 

experimental responses (A) were taken from a representative cell 
maintained in current-clamp at the indicated holding potential and 
stimulated either with current pulses or with mossy fiber bundle 
stimulation. Identical stimulation patterns were used for the model 
(B). (Ai) The traces on the left show the responses of a UBC 
maintained at — 80mV and injected with currents steps of different 
intensity (protocol below the voltage traces). Depolarizing current steps 
of adequate intensity generate a low-threshold spil<e (LTS) surmounted by 
a burst of sodium spikes. At higher intensity, the burst is followed by a 
tonic action potential discharge. With hyperpolarizing current steps, the 



UBC response shows a sag at the beginning of the step and a rebound 
burst at the end of the step. The trace on the right shows the responses 
of the same UBC maintained at — 68.2mV to a depolarizing current step. 
In this case the LTS and the associated burst almost disappear. (Aj) The 
traces on the /eft show the LOR evoked in the UBC [same as in (Ai)] 
by mossy fiber bundle stimulation with pulse-trains of different length. 
The traces on the right show the LOR evoked by mossy fiber bundle 
stimulation with pulse-trains of different intensity (in V). (Bi,B2) The 
traces show the responses of the UBC model corresponding to those 
generated in the experimental recordings shown in (Ai,A2). Note that 
the intensity of stimulation is in arbitrary units (low, medium, and high). 



(1) UBCs were silent at rest, with a resting potential of —66.7 ± 
1.8 mV (m = 81) (cf. -67.6 ± 0.9 mV in Diana et al, 
2007; -67.5 ± 1.4 mV in Locatelli et al, 2013). 

(2) During hyperpolarizing current pulses (from — SOmV), 
UBCs generated a sagging inward rectification followed 



by a rebound depolarization after the end of the pulse 
(Figure 2A1, left). 
(3) During depolarizing current pulses (from — SOmV), UBCs 
generated a spike burst, which, by raising current intensi- 
ties, was followed by a tonic discharge (Figure 2A1, left). 
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FIGURE 3 I Membrane mechanisms during repetitive firing induced by 
current steps. (A) Time course of ionic currents, calcium concentration, 
and membrane conductance of the model caused by a burst of action 
potential elicited by a depolarizing current step from two different holding 
potentials: left. +16 pA from -80 mV; right, +3 pA from -67.5 mV. The inset 
shows a detail of the persistent sodium current component. (B) Calcium 
currents measured separately in the soma and brush compartments of the 
model. (C) Time course of ionic currents taking part to rebound burst 
generation in the model following — 28pA current injection from — 80mV. 
Names and symbols as in Methods. 



The spike burst was composed of a LTS surmounted by a 
burst of high-frequency fast spikes. An AHP appeared at the 
end of the spike discharge. 
(4) When the resting membrane potential was artificially depo- 
larized (between —65 to — 70mV), the LTS almost disap- 
peared and the spike burst became weaker (Figure 2A1, 
Right). 

These electroresponsive properties were common to those 
reported in the previous experimental investigations (Diana et al, 
2007; Birnstiel et al, 2009; Locatelli et al, 2013), in which the 
implication of specific ionic channels was examined (see below). 
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FIGURE 4 I IVIembrane and cytoplasmic mechanisms during LOR 
generation. Time course of membrane and cytoplasmic mechanisms 
activated by a synaptic stimulus train (10 pulses at 100 Hz, arrows and 
vertical strip) delivered to the UBC model at the holding potential of 
—80 mV. Alpha is the forward rate function graded with synaptic inputs, 
which activates the intracellular cascade by converting the receptor X into 
X*. The amount X' determines the conversion of the intracellular factor Y 
into y* at rate Delta. Y* causes changes in H and TRP conductances and 
currents. The total current flux {1^ + Itrp) triggered by the cascade 
mechanism depolarizes the model generating the LOR. 



Following mossy fiber bundle stimulation in the presence 
of GABA and glutamate receptor antagonists, UBCs showed a 
LOR composed of a depolarizing ramp followed by a LTS sur- 
mounted by a high-frequency spike burst followed by a prolonged 
tonic discharge (Figure 2B1). The LOR corresponded to the same 
phenomenon reported by Locatelli et al. (2013), who suggested 
that the LOR depended on a metabotropic mechanism involving 
cAMP increase in the cytoplasm, leading to the H- and TRP- 
channel modulation. Notably, the delay, duration and frequency 
of the LOR were apparently related to the intensity and number 
of pulses used for mossy fiber bundle stimulation (Figure 2B1). 
Higher intensity or number of pulses caused shorter delay, higher 
frequency, longer duration of the LOR (quantitative data are given 
below in Figures 6, 7). 

The UBC model was parameterized against responses to step- 
current injection (Figure 2B1). The parameters related to spikes 
elicited by depolarizing current steps (e.g., using 16 pA from a 
holding potentials of — 80mV) included spike delay (30.3 ms), 
spike overshoot (36.7 mV), instantaneous frequency (176 Hz), 
steady state frequency (126 Hz), and frequency ratio (0.715). 
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FIGURE 5 I The effect of ionic current blocl<age on LOR generation. 

Simulations were made to mimic the pharmacological effect of ionic channel 
blockers on LOR generation. The LOR was elicited by a synaptic stimulus 
train (10 pulses at 100 Hz, arrows) delivered to the UBC model from the 



holding potential of —80 mV. The traces show the effect of 1^, Itrp, IcsHVA 
and combined Icshva + IcaLVA block. For comparison, the effect of In and Itrp 
block on membrane potential transients generated by hyperpolarizing current 
step injection is shown in the inset. 



The parameters related to hyperpolarizing current steps included 
sag amplitude (e.g., — 17.3 mV with a —20 pA step from —80 mV) 
and current threshold for rebound bursting (e.g., —20 pA 
from — 80mV). The match of model to the experimental data 
was statistically significant (Student's f-test, p < 0.01) for all the 
parameters. This model, once endowed with coupling mech- 
anisms regulating TRP and H channels, predicted the LOR 
following synaptic stimulation (Figure 2B2). 

It should be noted that what has been developed in this paper 
is a canonical model matching the properties of an ideal UBC. 
Variants in terms of burst length, frequency ratio, or rebound 
depolarization observed in some UBCs could all be obtained 
through slight changes of maximum current densities for specific 
ionic channels (data not shown). 

UBC ELECTROGENIC MECHANIMS IN THE MODEL 

In the model, as in real UBCs, a small depolarizing current 
(8-16 pA) from a negative potential (— 80mV) was sufficient to 
trigger the LTS and a short spike burst, while a larger depolariz- 
ing current (>16pA) was required to drive the model into tonic 
discharge (Figures 2B1, 3A; cf Locatelli et al, 2013). However, 
when depolarizing step-current injection was applied from rest- 
ing potential (-67.5 mV; Diana et al, 2007; Locatelh et al, 2013; 
Figure 2B1, Left), the response was characterized by a tonic dis- 
charge almost lacking the initial spike bust (Figure 2B1, Right). 
Even by using a net current injection stronger than that used 
from — 80mV, the instantaneous spike frequency was apparently 
lower. While a classical set of currents — including /jva, Ikv, ^ka, 
IcaHVAy and Ikcu — ^was sufficient to determine the spike shape 
as well as spike frequency and its variations (D'Angelo et al.. 



2001), specific ionic channels were required to determine these 
characteristic properties of UBC electroresponsiveness. 

The high responsiveness of the UBC model to small current 
injection from negative holding potentials was due not only to the 
relatively high input resistance but also to activation of the CaLVA 
current. CaLVA was the most prominent current in this functional 
regimen and regeneratively amplified the effect of current injec- 
tion causing the LTS. CaLVA currents have a voltage-dependent 
inactivation, reflected into the marked weakening of bursting 
elicited from more depolarizes membrane potentials (Diana et al., 
2007; Birnstiel et al, 2009; Locatelh et al, 2013). The density 
of CaLVA currents was derived from voltage-clamp recordings 
(Diana et al., 2007) and, together with a precise setting of input 
resistance, provided a major constraint for calibrating the whole 
electroresponsive mechanism. Once the LTS was generated, the 
depolarization gated the other ionic channels so that various cur- 
rents contributed to generate the spike bursts and control its 
evolution. As in real UBCs, burst spike frequency could raise 
up to 200 Hz along with spike amplitude adaptation. The burst 
was probably enhanced by the resurgent Na current (Raman and 
Bean, 1997; Afshari et al, 2004; Dover et al, 2010) and the first 
spike delay was regulated by Ika- 

With depolarizing current injections <16pA, the burst was 
driven by the LTS and terminated after a time compatible 
with CaLVA channel inactivation. With depolarizing current 
injections >16pA, spike discharge protracted beyond the burst 
through the intervention of another set of ionic channels: the 
CaHVA current and the persistent Na current sustained a depo- 
larizing plateau maintaining spike discharge. The model pre- 
dicted that an appropriate balance between the CaHVA and a 
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FIGURE 6 I Model robustness. The plots illustrate the results of varying the 
model maximum ionic conductances between the minimum and maximum 
values allowing to obtain existing granule cell responses. The target 
parameters were spike delay, instantaneous frequency, steady-state 
frequency, frequency ratio (Instantaneous/steady-state frequency), spike 
overshoot, and LOR duration. The polar plot reports the variations observed 
in each of these parameter by changing the maximum conductance of each 



slow K current, Kslow, was required to regulate the intensity 
and duration of this protracted spike discharge. Finally, dur- 
ing depolarization, the H current was deactivated and therefore 
reduced the total input conductance favoring depolarization, but 



ionic channel in turn. The histograms compare model to experimental results: 
stars represent the experimental data, the red dot represents the response 
of the canonical UBC model, color bars represents the variation of each 
parameter while changing the maximum ionic conductance between the 
extremes of the variation range. (A) Model robustness in response to step 
current injection (-1-16 pA for 800 ms from — SOmV). (B) Model robustness of 
the LOR. (C) Combined model robustness for step and LOR responses. 



its contribution to the total UBC current was limited by the 
decreased driving force. 

As in real UBCs, when injected with hyperpolarizing current 
steps, the model generated sagging inward rectification followed 
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FIGURE 7 I The LOR response space. Tlie LOR response space was 
reconstructed by systematically varying the number of pulses of mossy 
fiber stimuli. Experimental data were reported on the LOR response space 
for comparison. (A) The three panels show the response space for spike 
delay, burst duration, and instantaneous frequency of the LOR. The 
symbols represent experimental data sets in different UBC recordings, in 
which the number of stimulus pulses was varied. The intensity in these 
plots is the dependent variable obtained at the intersection between the 



number of mossy fiber stimuli and the parameter values measured 
experimentally and reported on the response space for that specific 
parameter (B) Ensemble intensity of mossy fiber stimuli obtained for each 
cell as the average of all values derived for the three plots in A (Mean ± 
MSE, n = 6 to 9 depending on the number of data points/plot). Dotted 
lines are linear fittings to the points of each data set. (C) Normalized SD of 
the intensity values calculated for each data set reported in (B) vs. the 
slope of the same data set. 



by rebound excitation after the end of the pulse. The rebound 
could be consisted of a depolarizing ramp that could prime a 
LTS and a spike burst when the hyperpolarizing was more intense 
and protracted (Figure 3C). The rebound was driven by H cur- 
rent deactivation and amplified by CaLVA current activation (see 
Figure 6; Russo et al, 2007; Locatelli et al., 2013). The H current 
was derived from voltage-clamp recordings (Locatelli et al, 2013) 
and provided a critical constraint for model calibration. 

MODELING THE LATE-ONSET RESPONSE 

By using the model, we tested the hypothesis that synaptic mod- 
ulation of H- and TRP-currents through an intracellular cascade 
was necessary and sufficient to generate the LOR (Locatelli et al, 
2013). Therefore, if the UBC model was correctly constructed, it 
should be able to generate the LOR by coupling synaptic trans- 
mission to H- and TRP-current gating through an appropriate 
second messenger pathway. The second messenger pathway was 
implemented as a generic reaction scheme, in which (1) mossy 
fiber bundle stimulation activated the extra-synaptic receptor X 
to X* , (2) X* increase above 60% of its maximum value initiated 



conversion of a second messenger 7 to Y*, and (3) Y* shifted 
the voltage-dependence of In and the maximum conductance 
of Itrp (see Methods; DiFrancesco and Mangoni, 1994; Locatelli 
et al., 2013; Figure 4). The increased Y* concentration eventually 
caused a proportionate increase in TRP maximum conductance 
and a proportionate shift of Ih activation toward positive mem- 
brane potentials. TRP depolarized the membrane sufficiently 
to generate the LOR, while Ih transiently unbalanced the sub- 
threshold current generating a slow depolarizing ramp tuning the 
delay of LOR initiation (Figure 4). It is worth noting that, after 
setting model intrinsic electroresponsiveness, no further changes 
was required to obtain the LOR, which is therefore a true model 
prediction. 

QUANTITATIVE COMPARISON OF THE MODEL WITH EXPERIMENTAL 
DATA 

The model could reproduce the fundamental properties of the 
LOR, which are observed in WCR recordings (this work and 
Locatelh et al., 2013). (i) The LOR was dependent on the inten- 
sity and number of synaptic pulses (cf Figure 2). (ii) The TRP 
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FIGURE 8 I Correlations of LOR parameters. Variation in LOR parameters 
vs. second messenger concentration. (A) The relative change of first-spike 
delay, burst duration and instantaneous frequency of the LOR were plotted 
against second messenger concentration. The first-spike delay decreases, 



while the burst duration and instantaneous frequency of the LOR increase 
with the cAMP level. (B) The relative change of different ionic currents with 
respect to second messenger concentration. Note the specific behaviors of 
H, TRR CaLVA, and KA currents. 



current generated the major depolarizing drive for LOR genera- 
tion, while the H-current regulated LOR delay by controlling the 
slope of the depolarizing ramp. Blocking H- and TRP-currents 
together abolished the LOR (Figures 4, 5). (iii) The Ca^+ currents 
enhanced LOR spike frequency, with a more evident contribution 
given by LVA than HVA currents. However, Ca^+ currents did not 
regulate the depolarizing ramp nor they were required to gener- 
ate the LOR itself (Figures), (iv) The LOR was poorly affected 
by holding potential and could be generated either from —80 
or —67.5 mV (data not shown). The reliability of model predic- 
tions about UBC electroresponsiveness and LOR generation were 
assessed through a quantitative comparison with experimental 
data. 

Ih provided a critical constraint for the model and the effect 
of pharmacological block by ZD7288 and Cs^ has been reported 
(Diana et al, 2007; Russo et al, 2007; LocateUi et al, 2013). Ih 
switch-off in the model caused hyperpolarization of resting mem- 
brane potential (Figure 5) and eliminated the sag during hyper- 
polarizing current steps (Figure 5, inset). The resting membrane 
potential shift (— 4mV) and differences between the maximum 
voltage deflection and the voltage attained at the end of the cur- 
rent step (from 3.96 to 0.13 mV) and the subsequent rebound 
depolarization (from 2.54 to 0.15 mV) were consistent with pre- 
vious experimental data (Locatelli et al., 2013). Finally, the LOR 
delay increased by 12.8% and the slope of the depolarizing ramp 
and its duration decreased by 12.2 and 24.9%, respectively, again 
reflecting experimental data (Locatelli et al., 2013). Therefore, 
in the model the Ih regulated the delay of the late-onset burst 
by controlling the depolarizing ramp similar as in whole-cell 
recordings. 

The TRP channels provided a background leakage current, 
which increased during the LOR. Itrp switch-off in the model 
prevented full LOR generation leaving only a subthreshold depo- 
larization driven by the H current (Figure 5). The slope of 
the depolarizing ramp decreased by 69%, consistent with 49% 



decrease measured experimentally following channel blockage 
with SKF96365 (Locatelli et al, 2013). 

Pharmacological blockage of CaLVA channels (with 
mibefradil) and CaHVA channels (with nimodipine) was 
reported to modify spike firing in the LOR but not to prevent 
LOR generation (Locatelli et al., 2013). In the model, CaLVA 
switch-off did not prevent the LOR but reduced the number 
of spikes by 20% and the instantaneous frequency by 38%. 
Combined CaHVA and CaLVA currents switch-off did not 
prevent the LOR but reduced the number of spikes by 80% and 
the firing frequency by 88% (Figure 5). These effects in the model 
are comparable to those of pharmacological blockage reported 
experimentally. 

ROBUSTNESS OF MODEL PREDICTIONS 

As explained above (see also Methods), reference values were 
obtained from UBC voltage-clamp data on the H- LVA- HVA- and 
TRP-currents and were then used to pre-set the corresponding 
maximum conductance values in the model. The other maximum 
conductance values were set respecting proportions evaluated 
in previous models [for example the Gmax (/Na)/Gmax i^Kv) 
ratio, the G^ax {lKCa)/G!n^{IcaHVA] ratio, (Gabbiani et al., 1994; 
D'Angelo et al., 2001; Solinas et al, 2007b). Nonetheless, the 
maximum ionic conductances in the model were actually free 
parameters, as they were tuned to match UBC electrorespon- 
siveness. A robustness analysis was then performed in order to 
estimate the confidence intervals of the maximum ionic con- 
ductances still allowing to generate existing UBCs. The range of 
variations allowed for each maximum conductance were used to 
maintain the model within the experimental values of first spike 
delay, instantaneous frequency, steady-state frequency, frequency 
ratio and spike overshoot (Figure 6). Finally the model predicted, 
through the cytoplasmic coupling mechanism, the experimental 
values measured in the LOR for the parameters indicated above 
as well as for LOR duration. 
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The experimental parameters measured in the UBCs covered 
in this study were used to set the range limits for parameter varia- 
tion in the model. Intrinsic excitability was characterized by spike 
delay (31.1-93.4 ms, « = 12), spike overshoot (15.0-45.1 mV, n = 
12), instantaneous frequency (95.9-202.7 Hz, n= 12), steady- 
state frequency (41.2-141.8 Hz, n = 12), and frequency ratio 
(0.61-1.3, n = 12). The LOR was characterized by spike delay 
(98-1000 ms, M = 18), spike overshoot (10.2-44.5 mV, « = 18), 
instantaneous frequency (92-209 Hz, « = 18), steady-state fre- 
quency (19.5-60 Hz, « = 18), frequency ratio (0.1-0.6, n = 18), 
and burst duration (100-2500 ms, n = 18). Simulation results 
were accepted as possible UBC responses when falling within 
these limits. 

Both concerning responses to current injection (Figure 6A) 
and the LOR (Figure 6B), the modulation of maximum ionic 
conductances led within boundaries spanning < ±20%. By com- 
bining the two sets of boundaries obtained for current injec- 
tion and LOR gave the compound boundaries (Figure 6C). This 
robustness analysis showed that the canonical model was well bal- 
anced and could use the same ionic channel settings to faithfully 
explain both electroresponsiveness and LOR. 

PREDICTION OF LATE-ONSET RESPONSE ACTIVATION BY DIFFERENT 
ACTIVITY PAHERNS 

The model allowed to predict burst delay, duration and frequency 
of LOR in relation to different combinations of stimulus dura- 
tion and intensity (Figure 7A). The response showed opposite 
changes for delay (which decreased with stimulus duration and 
intensity) compared to duration and frequency (which increased 
with stimulus duration and intensity). If the simulated response 
space was correct, the experimental data points obtained chang- 
ing stimulus duration should actually correspond to the same 
stimulus intensity. Actually, for all parameters, the experimental 
data points fell on nearly horizontal lines (Figure 7A) . A quantita- 
tive assessment of the matching between experiments and model 
was obtained by calculating the slope of the duration/intensity 
relationships (Figure 7B) and by assessing the variability of the 
mean predicted stimulus intensity (Figure 7C). The slope of 
the duration/intensity relationships was 0.0005 ± 0.002 and was 
statistically indistinguishable from 0 {p < 0.001, n = 7, f-test). 
The variability of the predicted stimulus intensity was <10%. 
Therefore, the model could effectively predict LOR properties 
while varying the stimulus pattern. 

PREDICTION OF MECHANISMS DETERMINING LOR PROPERTIES 

The explanation of why, in the LOR stimulus/response space, 
similar parameter values were determined by varying either stim- 
ulus duration or intensity (cf Figure 8A) had to be searched in 
their control over the common factor, Y*. In turn, Y* controlled 
membrane channel gating and therefore the LOR (Figure SB). 

The effect of Y* on LOR parameters (burst delay, duration 
and frequency) for all possible combinations of stimulus dura- 
tion and intensity used in Figure 7 are shown in Figure 8A. 
In particular, while increasing Y*, burst delay decreased almost 
exponentially, burst frequency increased almost exponentially, 
and duration increased almost linearly (Figure 8A; the pres- 
ence of multiple points at each Y* value reflected the different 



Y* kinetics for different combinations of stimulus intensity and 
duration). 

A mechanistic explanation of the LOR properties was obtained 
by correlating Y* with the underlying ionic currents (Figure SB). 
It turned out that the Y* dependence of I^a, Ikv, IksIow, IcaHVA^ 
iKCa, and IjRp was similar to that of LOR duration. The Y* depen- 
dence of Ih and Ika was similar to that of LOR frequency and 
opposite to that of LOR delay. The Y* dependence of IcaLVA was 
inverted compared to any others. This result implies a complex 
interaction of all currents in regulating LOR properties and a 
primary role of In in controlling LOR delay. 

DISCUSSION 

This paper shows that regulation of H- and TRP-currents by an 
intracellular factor can generate the LOR, a slow depolarization 
driven by synaptic activity recently observed in cerebellar UBCs. 
The LOR mechanism was explained by a model accounting for 
UBC electroresponsiveness and including cytoplasmic coupling 
to H- and TRP channels. Compared to glutamatergic synaptic 
responses, the LOR generated a spike burst with longer delay, 
which was modified by the pattern of mossy fiber discharge. The 
LOR, together with specialized mechanisms of intrinsic electrore- 
sponsiveness and synaptic transmission, endows UBCs with the 
ability to generate spike bursts with variable delay in the cerebellar 
granular layer. 

THE UBC CANONICAL MODEL: PROPERTIES AND LIMITS 

The present model provides a biologically realistic reconstruc- 
tion of UBC electroresponsiveness (Diana et al., 2007; Russo 
et al., 2007; Birnstiel et al, 2009; Locatelli et al, 2013) based 
on the representation of neuronal geometry, passive proper- 
ties and on the subsequent insertion of appropriate ionic and 
intracellular mechanisms. The H, TRP, and Ca^+ conductances 
provided the critical constraints for UBC modeling and had den- 
sity, gating and kinetics derived from UBC recordings (Diana 
et al, 2007; Russo et al, 2007; Birnstiel et al, 2009; LocateUi 
et al, 2013). These core mechanisms were associated with a 
spike generation mechanism sustained by Na^" and conduc- 
tances. The CaLVA and CaHVA channels were distributed across 
soma, dendrite and brush in order to generate calcium cur- 
rents proportional to those determined experimentally (Diana 
et al, 2007; Birnstiel et al., 2009), while the Na channels were 
localized in soma, initial segment and axon. The localization 
of ionic channels had an impact on model electrogenesis. On 
one side, the CaLVA current was larger in the brush than 
soma, while the CaHVA current was equally distributed, as indi- 
cated by calcium current measurements (Diana et al, 2007; cf 
Figure 3C). The high density of CaLVA channel in the brush, 
coupled with limited current flow to the soma through the 
dendrite, enhanced the local auto-regenerative LTS depolariza- 
tion. This would become particularly effective when synaptic 
currents are injected into the brush. On the other side, the 
Na+ current was concentrated in the soma, initial segment and 
axon, conferring higher excitability and normalizing spike and 
burst discharge parameters (not shown). The precise localiza- 
tion of Na+ channels in UBCs remains to be determined by 
immunohistochemistry. 
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The resulting model was a canonical (or prototypical) rep- 
resentation of the UBC, such that ±20% variations of ionic 
channel maximum conductances still allowed model responses 
(first spike delay, initial firing rate, steady-stated firing rate, fre- 
quency ratio, spike overshoot) to remain within the range of 
parameters recorded experimentally. This variation range is nor- 
mally observed while measuring ionic currents in neurons (Bhalla 
and Bower, 1993) and has also been reported in UBCs for H- 
, CaLVA, and CaHVA currents (Diana et al, 2007; Birnstiel 
et al., 2009; Locatelli et al., 2013). Most importantly, the model 
allowed to predict the LOR without any further modifications, 
provided that the cytoplasmic coupling mechanism was properly 
implemented. 

Beside its effectiveness, some aspects of the UBC model were 
admittedly simplified, (i) The H-current kinetics observed in 
UBC recordings (Locatelli et al, 2013) suggest that H-current 
gating could reflect differential expression and assembly of sub- 
units. This could contribute to tune UBCs toward specific LOR 
kinetics, (ii) The Na current was implemented using a 13-state 
stochastic model for Navl.6 channels of granule cells generating 
transient, persistent and resurgent current components (Raman 
and Bean, 2001; Magistretti et al, 2006; Dover et al, 2010). 
While these three components are present in all cerebellar neu- 
rons (Afshari et al., 2004; Khaliq and Raman, 2006), the resurgent 
current may be proportionally higher in UBCs than granule cells 
(Afshari et al, 2004). This could further promote the initial spike 
burst primed by the LTS. (iii) The regulation of duration and fre- 
quency of tonic discharge following the LTS required that a slow 
K current, Kslow(D'Angelo et al., 2001), was added to the model. 
Moreover, the regulation of first-spike delay required that a fast- 
activating K current, KA(D'Angelo et al, 2001), was added to the 
model. These current have not been reported in UBCs yet and are 
therefore model predictions, which remain to be proved experi- 
mentally, (iv) Intracellular calcium regulation was modeled using 
a simplified scheme useful to couple calcium entering through 
CaHVA channels to BK-type channel gating and the effect of more 
complex schemes remains to be explored. 

In a further evolution, the model may benefit of a detailed 
morphological reconstruction of UBC brush and axon along with 
proteomic and channelomic analysis directly defining the molec- 
ular properties of ionic channel subtypes and of the cytoplasmic 
cascade coupling membrane receptors to ionic channels. It will 
also be of interest to determine under which circumstances rou- 
tines based on genetic algorithms will be able to solve the complex 
optimization problem of the UBC model (Druckmann et al., 
2007, 2008, 2011,2013). 

THE THREE CORE MECHANISMS OF BURST GENERATION IN UBCs 

A quite relevant property of the UBC model, as well as of real 
UBCs, was that of generating a variety of burst responses, which 
were controlled by the intervention of H-, TRP-, and Ca^^"- cur- 
rents. The model showed that the interaction of H- TRP-, and 
Ca^"*"- currents emerged in three conditions. 

(1) Bursts during depolarization from negative memhrane poten- 
tial. In response to depolarizing current injection, the IcaLVA 
was critical to promote spike generation. 



(2) Bursts at the end of a hyperpolarization. Following a hyper- 
polarization, Ih caused a rebound excitation boosting the 
voltage-dependent activation of IcaLVA and promoting gen- 
eration of rebound bursting. 

(3) Bursts during the LOR. During the LOR, the ramp was pro- 
moted by Itrp and regulated by Ih. IcuLVA then amplified the 
depolarization intensifying the subsequent spike burst. 

These mechanisms, especially the first two, proved markedly 
voltage-dependent due to the specific voltage-dependent inacti- 
vation of CaLVA channels. It is therefore predicted that alternat- 
ing cycles of depolarization and hyperpolarization involving the 
granular layer inhibitory circuit (Vos et al., 1999) could modify 
the ability of the UBCs to generate spike bursts substantially. 

LOR GENERATION IN THE MODEL 

Rather than reproducing the intracellular biochemical path- 
way (which remains largely undetermined) we have developed 
a generic coupling mechanism including receptors activated by 
synaptic activity and leading to the production of a second mes- 
senger (y*), possibly corresponding to ATP/cAMP conversion. 
The level of this second messenger caused the current influx 
(~-16pA) generating the depolarizing ramp and the LOR, with 
a delay depending on the second messenger concentration and 
ranging up to several hundreds of milliseconds. 

Thus, the model provided substantial support to the existence 
of an intracellular mechanism coupling membrane receptor acti- 
vation to H- and TRP-channels and causing the LOR in UBCs. 
Modeling suggests that In was responsible for controlling LOR 
delay and frequency, that IcaLVA and Ika were driven into the 
control of LOR delay and frequency through their time- and 
voltage-dependent inactivation, that I^a, IcaHVA and Ij-rp cur- 
rents drove the depolarization during the burst and determined 
burst duration and activation of Ikv> IksIow> and lKCa> which 
actually regulated the frequency and terminated the burst. 

THE LOR AND DELAY-LINES IN THE CEREBELLUM 

Through the LOR, UBCs can translate the intensity of the input 
(coded as the number of active fibers and duration of their dis- 
charge) into output bursts with different delay, duration and 
frequency. This implements a time-code that could reverberate 
through the network along the chains made by UBCs with other 
UBCs and granule cells (Nunzi and Mugnaini, 2000; Nunzi et al., 
2001). The LOR could be modulated by local network activ- 
ity causing various patterns of UBC excitation and inhibition 
(Ruigrok et al., 201 1; Rousseau et al., 2012). Recently, a frequency- 
dependent regulation of AMPA receptor-mediated after responses 
has also been reported (van Dorp and De Zeeuw, 2014) and 
has been suggested to provide a system to generate protracted 
responses in UBCs. The implementation of these AMPA receptor- 
dependent mechanisms into the UBC model would allow to test 
its actual impact on UBC responses and its relationship with 
the LOR. 

It has been proposed that UBCs are essential for shifting and 
converting the phase of mossy fiber activity that relays infor- 
mation from the vestibular apparatus, eyes, or neck, and that 
their characteristic cellular properties are particularly relevant 
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for controlling and consolidating motor learning in paradigms 
such as VOR phase reversal (Gao et al., 2012). Actually, the LOR 
could help implementing the "velocity storage" system (Dai et al., 
2007), which allows the vestibulo-cerebellum to transform the 
head velocity signal into commands for ocular motor neurons 
controlling the slow phase of nystagmus (Goldberg et al., 1984; 
Okada et al., 1999). Patterns of activity compatible with the inter- 
vention of the LOR in UBCs have been reported in the electric 
organ (a cerebellar-like structure) of mormyrid fishes during the 
prediction of the sensory consequences of motor acts (Kennedy 
et al., 2014). By involving intracellular cascades, the LOR could 
be modified by receptor systems — involving e.g., acetylcholine, 
serotonin, and noradrenaline — known to modulate H and TRP 
channels, thus correlating cerebellar computation with the brain 
functional state. 

CONCLUSIONS 

The UBCs, despite their relatively recent discovery, are revealing 
specific cellular mechanisms generating a rich set of electrore- 
sponsive properties. By these means, UBCs generate different 
types of burst responses including the LOR. The LOR endows the 
granular layer with a long-sought slow process extending the tim- 
ing capabilities of the cerebellar network (Dean and Porrill, 2008). 
The model predicts that the H-current and the TRP-current are 
critical to generate the LOR through membrane receptors and 
cytoplasmic cascades that remain to be elucidated. The model also 
predicts the specific roles for the different types of Ca^+ and Na 
currents observed in UBCs and predicts the existence of fast and 
slow-K currents (of the category oflKA and IksIow)' which remain 
to be demonstrated experimentally. The UBC model provides a 
useful new tool for investigating granular layer spatio-temporal 
dynamics and for implementing large-scale network models of 
the cerebellum. 
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